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We study the quantum interference effect for the single ballistic Aharonov-Bohm billiard in 
the presence of a weak magnetic field B. The diagonal part of the wave-number averaged reflection 
coefficient 81Zd is calculated by use of semi-classical scattering theory. In addition to the appearance 
of "h/2e oscillation" that are caused by interference between time-reversed coherent backscattering 
classical trajectories, B in the conducting region leads to negative magnetoresistance and dampening 
of the h/2e oscillation amplitude. The B dependence of the results reflects the underlying classical 
(chaotic and regular) dynamics. 
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I. INTRODUCTION 



Recently, the interplay of chaos and quantum interference in ballistic quantum dots has intrigued experimentalists 
and theorists alike. The quantum interference effects, e.g., ballistic weak localization (BWL) |l| and ballistic conduc- 
tance fluctuation Q in such structures depend on whether the underlying classical dynamics is regular or chaotic. 
Therefore, these effects are interesting from the viewpoint of the theory of quantum chaos || . 

More recently, we predicted that h/2e oscillation of magnetoconductance, analogous to Altshuler-Aronov-Spivak 
(AAS) effect [Q in disordered systems ||, should be observable in a single ballistic Aharonov-Bohm ring (hereafter 
called AB billiard) with magnetic flux penetrating only through the hollow This phenomenon of conductance 
oscillation is caused by the interference between a pair of time-reversed coherent back-scattering classical paths that 
wind the center obstacle in the billiard. We calculated the diagonal part of the BWL correction to the wave-number 
averaged reflection coefficient by use of semiclassical scattering (SCS) theory . Our analysis |j| yielded for the 

chaotic AB billiard 



5Tl D ($) ~ Vexp(-<fo)cos (4tto— ), (1) 

n— 1 v ' 

where &g = h/e is the magnetic flux quantum and <f> is the magnetic flux that penetrates the hollow. In eq. (^) 
S = y / '2To'j / a, where a, To, and 7 are system-dependent constants and correspond to the variance of the winding 
number distribution ||, the dwelling time for the shortest classical orbit and the escape rate @J§], respectively. In 
this case, the oscillation amplitude decays exponentially with increasing rank of higher harmonics n, so that the main 
contribution to the conductance oscillation comes from the n=l component that oscillates with the period of h[2e. 
By contrast, for regular and mixed (Kolmogorov-Arnold-Moser tori and chaotic sea) AB billiard, we obtained 



6K„(Z) ~£F^--,z+-;-^-jcos [ 4 ™t ) > (2) 

n— 1 x / \ / 

where F and [3 are the hypergeometric function of confluent type and the exponent of dwelling time distribution 
N(T) ~ T~ z 0,^), respectively. In eq. (|J) the oscillation amplitude decays algebraically for large n, and therefore 
the higher-harmonics components give noticeable contribution to magnetoconductance oscillations. These discoveries 
indicate that the h/2e AAS oscillation occurs in both ballistic and diffusive systems forming the AB geometry and the 
behavior of higher harmonics components reflects a difference between chaotic and non- chaotic classical dynamics. 

In real experiments, however, the magnetic field would be applied to the entire region (both the hollow and annulus) 
in the billiard. Thus it is indispensable to apply the SCS theory to this case, in order to see the comparison between 
the experimental data and theoretical prediction. In this situation, we shall envisage h/2e oscillation together with 
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the negative magnetoresistance and dampening of the h/2e oscillation amplitude with increasing magnetic field. In 
this paper, we shall focus our attention on two-dimensional ballistic AB billiards (e.g., the insets in Fig. 1) with the 
magnetic flux penetrating through the entire region and calculate reflection amplitude by use of SCS theory. 



II. SEMICLASSICAL THEORY 







Following Baranger, Jalabert and Stone's arguments (JJffl, we start with a quantum-mechanical reflection amplitude 



S n , m ~ ihyjv n v m dy dy'ipl(y')ilj m (y)G(y' ,y,E F ), 



(3) 



where v m (v n ) and V'm(V'n) are the longitudinal velocity and transverse wave function for the mode m(n), respectively. 
G is the retarded Green's function. To approximate r n . m we replace G by its semiclassical Feynman path-integral 
expression jlll , 



G sc (y',y,E) 



2tt 



{2-KiTif/ 2 



^2 \ZoTexp 



s(y,y') 



^S s {y',y,E) ~ i^n s 



(4) 



where S s is the action integral along classical path s, D s = (v F cos9')~ 1 \(dO/dy') y \ , 9 (9') is the incoming (outgoing) 
angle, and /i s is the Maslov index. Assuming hard walls in the leads, we substitute eq. (Q) into eq. (||) and carry out 
the double integrals by a stationary-phase approximation. Thus we obtain 



2W 



sgn(n)sgn(m) y D s exp 



s(n,m) 



Z ~ 7T 

-S s (n,m;E) - i-p, s 



(5) 



where W is the width of the hard-wall leads and m = ±m. The summation is over trajectories between the cross 
sections at x and x' with angle sin 6 = nn/kW. In eq. (||), S s (n,rh; E) = S s (y' Q ,yo; E) + hir(fhyo — ny' )/W, 
D s = (m e VF cos9')~ 1 \{dy/d9')o\ and jl s = fi s + u (— (d9/dy)' y ) + u(—(d9' /dy')g) , respectively, where u is the 



Heaviside step function. The Kronecker delta term in eq 
of zero length |Q . Within the diagonal approximation 
probability R c i, viz., 



is exactly canceled by the contributions of paths 
the quantum correction SR to the classical reflection 



R 



N M 

y ] \ r n,m\ 
n,m— 1 



Rci + 5R 



with the mode number Nm, is given by 

SRp = ^ ^ \l A S A U exp ik(L s -L t 



(6) 



(7) 



where s and u label the classical trajectories. In eq. (R), L s — S s /k% , v SiU — (ji u — \i s ) /2 , and A s = (hk/W) D s 
. The wave-number averaging of SRd over all k, denoted as SIZd , eliminates all paths except those that satisfy 
L s = L u in eq. (^). In the absence of spatial symmetry, L s = L u holds if u is the time reversal of s. A weak 
magnetic field does not change the classical trajectories appreciably but does change the phase difference between the 
time-reversed trajectories by (S s — S u )/h = 20 s B/&o, where 6 S = 2ir j s A ■ d£/B is the effective area almost enclosed 
by the classical path. To evaluate the summation over s and n, we shall reorder the backscattering classical paths 
according to the increasing effective area. Therefore, we obtain 



sn D {B) 



dON(e) exp 



20B\ 



(8) 



where N(Q) is the distribution of 0. The phenomelogical statistical theory leading to a distribution of the enclosed 
area N(Q) for chaotic AB billiards is given as follows. There exist two kind of classical paths and N(Q) is the sum 
of the distribution of unwinding trajectories, viz., Nq(Q) and that of 0) winding trajectories. This is due to 
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all classical trajectories winding a center obstacle n times until they exit, except for very short backscattered paths 
(n=0 component). iVo(0) is essentially the same as that of an ordinary chaotic billiard (e.g., stadium), obeying a 
monotonic exponential law |@j9j, i.e., iV (9) ~ exp(— £ c \ |0|) , where e c i is the inverse of the average area enclosed 
by classical trajectories. Therefore, the full distribution of the enclosed area is given by 



N(e)~N Q {<d)+ N(e,n)P(n), 



(9) 



where P(n) and N(Q,n) are the distribution of the winding number n and that of the enclosed area for a given 
n, respectively. Owing to the ergodic properties of fully chaotic systems, N(Q,n) is assumed to obey a Gaussian 
distribution in which the variance of area is proportional to n , i.e., 



N(Q, n) 



1 



exp 



2/3 \n\ 



On the other hand, exploiting Berry and Keating's argument g, P(n) is given by 



P{n) 



/ dTP(n, T)N(T) ~ exp ( 
Jo 



-S\n\), 



where S = y/2T j/a. In eq.(|llj), N(T) ~ exp(- 7 T) and 

P(n,T)-- 



T 



2naT 



exp 



n 2 T, 



2aT J 



(10) 



(11) 



(12) 



are the exponential dwelling time distribution [[7|J9[| and the Gaussian distribution of winding numbers n for trajectories 
with a fixed T [pUq], respectively. With the use of eqs.(llG) and (O), we reach 



JV(0) - Ae^'l ! + 



1 



: exp 



2/3 \n\ 



-s\ 



(13) 



To examine the validity of expression (|T^), we directly calculated N(Q) for chaotic AB billiards (a single Sinai billiard 
|n|) by mean of classical numerical simulations. In the calculations, we inject 10 8 particles into the billiard at different 
initial conditions. N(Q) for the chaotic AB billiard has proved to be nicely fitted by eq. ( |l3| ) [see Fig. 1(a)]. (In 
this case, ©o/2-7r is approximately the average area between the outer square and inner circles.) As the size of center 
obstacle approaches zero, the Sinai geometry becomes square (regular billiard). In this case, we have confirmed that 
the oscillation structure disappears and JV(0) obeys a well-known power law |7]|J. Substituting eq. ( p^ ) into eq. (^), 
we finally obtain 



6K D (<5>) 



1 



/ 47T <I> 

\OqScI $0 



exp 

n=l 



- /3/47T $ 



cos 47m- 



(14) 



where $ = BQq/2it. The first term in eq. (|l4|), i.e., Ae^ 1 / |l + (2_B/e c ;$o) 2 | ; which contributes to negative 
magnetoresistance, agrees with Baranger, Jalabert and Stone's Lorentzian BWL correction 



5K D {B) 



1 + 



2B 



(15) 



for chaotic billiard m\M , where 1Z c i is the wave- number-averaged classical reflection probability. 
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FIG. 1. Semi-logarithmic plot of the effective area distributions in scattering from the (a) Sinai (chaotic) billiard (R/W = 3 
and L/W = 10) and (b) square AB (regular) billiard (LinjW = 6 and L ut/W = 10). The numerical simulation results 
(diamond) for (a) and (b) are well fitted by eqs. (13) and (16) (solid line), respectively. The insets show the schematic views 
of two types of AB billiards. 

For the second term in eq. (|T^) (n > components), the oscillation amplitude decays exponentially with increasing 
n. Therefore, the main contribution to the conductance oscillation comes from the n=l component, which oscillates 
with period h/2e. This behavior is consistent with our previous specific result [i.e., eq. ([!])] for the chaotic AB billiard 
in which magnetic flux penetrates only through the hollow [0] . In addition to this property, the oscillation amplitude 
damps exponentially with increasing magnetic field. 

On the other hand, for regular cases (the square AB billiard) the form of N(Q) has been estimated as 



AT(e)~n (|e|+A 1 p 2 + ]T (|n|+n 1 )^ 1 (|8-ne |+A 2 )- 



(16) 



from the numerical simulation [see Fig. 1(b)]. In the calculation wc injected 9 x 10 8 particles into the billiard. In cq. 
( p"6| ) no, Tlx, z\, Z2, Ai, and A 2 are also fitting parameters and 0q/27t is approximately the average area between the 



outer and inner squares in this case. This distribution leads to |14 



5K D (<f>) ~ n Ai(Ai, $) + 2Ai (A 2 , $) (n + m) 21 cos Unn 



where 



71=1 



Ai(A,$)= / dx(x + AY 
Jo 



cos I X 

6 $ 



(17) 



(18) 



Since Ai(Ai,$) is equal to the Fourier transform of a power-law function, one can expect a cusplike BWL peak 
near zero magnetic field [Q,^5[. In contrast to chaotic cases, the oscillation amplitude decays algebraically for n. 
Therefore, we can see that higher-harmonics components give a significant contribution to conductance oscillations. 
This is because the number of multiple-winding trajectories is much larger in regular billiards than in chaotic billiards. 



III. AAS OSCILLATION AND NEGATIVE MAGNETORESISTANCE 

In this section we shall discuss in more detail the difference of 61Zjj{&) between chaotic and regular AB billiards. 
In Fig. 2 we show S1Zd(&) for Sinai (chaotic) and square AB (regular) billiards. The values of the fitting parameters, 
determined by the classical simulation, are substituted into eqs. (|l4]) and ([l~7]). In order to see the marked difference 
of the <& dependence of S1Zd(®), we shall investigate the n = term in eqs. (|l4|) and fll7j), denoted as 61Znmr{$) 
and the n > term, denoted as S1Zaas(®), separately, i.e., 
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FIG. 2. Semi-classical BWL correction 81Zd to the reflection coefficient for chaotic AB (Sinai) billiard (solid) and regular 
square AB billiard (dotted) as a function of 4>(= QoB/2-k). SIZd is normalized to the value at $ = , i.e., the classical 
reflection probability TLd- 

sn D {<s>) = 5k NM r{$) + sn AAS {<p). (19) 

Figure 3(a) shows 61Znmr(&) which contributes to the negative magnetoresistance for two types of billiards. The 
shapes of 51Znmr,(&) in the vicinity of zero magnetic field are quite different between chaotic and regular billiards, 
i.e., a quadratic curve versus a linear line [see the inset in Fig. 3(a)]. For large $ (but with the cyclotron radius 
sufficiently large compared to the system dimension), 51Znmr(&) saturates in chaotic cases, but shows nosaturation 
in regular cases. Similarly, the STZaas ($) corresponding to the h/2e AAS-like oscillation part is indicated in Fig. 3(b). 
While the oscillation amplitude damps rapidly with increasing $ for the chaotic AB billiard, it damps gently for 
regular AB billiards. 

Therefore, on the basis of the above results, the qualitative difference of S1Zd(&) between chaotic and regular AB 
billiards is attributed to the different classical distribution of the effective areas. As the dimension of the center 
obstacle (e.g., R for a Sinai billiard and Li n for a square AB billiard) becomes zero, the oscillation structure of N(Q) 
is indistinct for two types of billiard, so that the h/2e conductance oscillation would disappear. 

To consolidate the above semiclassical prediction of the h/2e oscillation, we must compare eqs. (0) and (0) with 
quantum- mechanical calculations (e.g., a recursive Green's function method p6| ) and also check the influence of the 
off-diagonal contribution to 61Z(Q) [p|[7|JT^]. Moreover, it is desirable to confirm our prediction by having recourse to 
a random matrix approach for systems with broken time reversal symmetries |Tg| ] . Such investigations will be given 
elsewhere. 

IV. CONCLUSION 

In summary, we have derived the semiclassical formula for 8Hd{§) of single chaotic and regular AB billiards in 
which a weak B is applied to the entire region. We have shown that h/2e oscillations and negative magnetoresistance 
would appear concurrently in 51Zd{<&)'- as for h/2e conductance oscillations, we find the oscillation mainly with a 
fundamental period h/2e and rapid damping of the amplitude with increasing B for the chaotic billiard versus the 
large contribution of higher harmonic components and mild damping of the oscillation amplitude for a regular billiard 
; As for negative magnetoresistance, the Lorentzian peak and saturation for the chaotic billiard versus a cusplike 
structure and no saturation for a regular billiard are reproduced. Although Taylor et al. recently made an 

experimental study of the weak localization peak and self-similar structure of magnetoresistance in a semiconductor 
Sinai billiard, no detailed experimental result of h/2e conductance oscillations has yet been reported. We hope that 
these characteristics of quantum chaos in the quantum magnetotransport will be experimentally observed in ballistic 
quantum dots forming AB geometry. 
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FIG. 3. Magnetic flux dependence of two different components of BWL corrections: (a) STZnmr contributing to negative 
magnetoresistance and (b) STZaas contributing to the h/2e oscillation for chaotic (solid) and regular (dotted) AB billiards. 
STZnmr and STZaas are normalized to the value at $ = 0. The inset in (a) shows STZnmr in the vicinity of zero magnetic flux. 
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